Transport properties of granular metals at low temperatures 
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We investigate transport in a granular metallic system at large tunneling conductance between 
the grains, qt 2> 1. We show that at low temperatures, T < gvS, where 5 is the single mean energy 
level spacing in a grain, the coherent electron motion at large distances dominates the physics, 
contrary to the high temperature (T > grS) behavior where conductivity is controlled by the 
scales of the order of the grain size. The conductivity of one and two dimensional granular metals, 
in the low temperature regime, decays with decreasing temperature in the same manner as that 
in homogeneous disordered metals, indicating thus an insulating behavior. However, even in this 
temperature regime the granular structure remains important and there is an additional contribution 
to conductivity coming from short distances. Due to this contribution the metal-insulator transition 
in three dimensions occurs at the value of tunnel conductance g^ = (I/Gtt) \n{Ec/S), where Ec is 
the charging energy of an isolated grain, and not at the generally expected oc 1. Corrections 
to the density of states of granular metals due to the electron-electron interaction are calculated. 
Our results compare favorably with the logarithmic dependence of resistivity in the high- Tc cuprate 
superconductors indicating that these materials may have a granular structure. 

PACS numbers: 73.23Hk, 73.22Lp, 71.30.+h 



A great deal of research in the current mesoscopic 
physics focuses on understanding properties of granular 
metals (see The interest is motivated by the 

fact that while their properties are generic for a wealth of 
strongly correlated systems with disorder, granular met- 
als offer a unique experimentally accessible tunable sys- 
tem where both the interaction strength and degree of 
disorder can be controlled. 

The key phenomenon revealing the most of the under- 
lying physics is transport, where the effects of interac- 
tions play a crucial role. The processes of electron tun- 
neling from grain to grain that govern electron transfer, 
are accompanied by charging the grains involved after 
each electron hop to another grain. This may lead to a 
Coulomb blockade, and one justly expects this effect to be 
of the prime importance at least in the limit of weak cou- 
pling. It makes it thus clear, on a qualitative level, that 
it is the interplay between the the grain-to-grain coupling 
and the electron-electron Coulomb interaction that con- 
trols transport properties of granular metals; yet, despite 
the significant efforts expended, a quantitative theory of 
transport in metallic granular systems is still lacking. 

A step towards formulation such a theory was made 
recently in [ Q . It was shown that depending on the 
dimensionless tunneling conductance gr one observes ei- 
ther exponential-, at (?t ^ 1, or logarithmic, at S> 1 
temperature dependence of conductivity. The considera- 



tion in [I3II was based on the approach developed by Am- 
begaokar, Eckern and Schon (AES) ^ for tunnel junc- 
tions. This technique however, as shown in 0, applies 
only at temperatures T > qt^, where S is the mean en- 
ergy level spacing in a single grain, in this regime the 
electron coherence does not extend beyond the grain size. 
At low temperature region, T < qtS, the effects of the 
electron coherent motion at distances much exceeding the 
single grain size a must be included, thus this important 
regime is not described by the AES approach ^ . 

Although experimentally the low temperature regime 
is well within the experimental reach P, Q] ) it has never 
been addressed theoretically so far. The important ques- 
tion whether the system is a metal or becomes an in- 
sulator, in other words, whether the conductivity of the 
granular metals at large conductances remain finite in 
the limit of T ^ is still open. 

In this Letter we investigate the low-temperature con- 
ductivity of granular samples focusing on the case of large 
tunneling conductance between the grains, (Jt ^ 1- To 
this end we develop a technique that goes beyond the 
AES approach and includes effects of coherent electron 
motion at distances larger than the size of the grain. 
Without the Coulomb interaction the granular system 
would be a good metal in the limit, ^ 1, and our 
task is to include the charging effects into the theory. 
We find that at temperatures, T < grS properties of the 
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granular metal depend on the dimensionality of the array, 
and corrections to the conductivity and density of states 
due to the effects of Coulomb interaction are similar to 
those obtained in Ref. for a homogeneous metal. Thus 
at low temperatures the systems behaves essentially as a 
homogeneous metal contrasting the case of large temper- 
atures, T grS considered in Ref. 0]. 

This means that at large conductances the 3D system 
is a good metal. On the other hand, at ^ 1 a granular 
sample is in the insulating state. Therefore a 3D system 
should exhibit a metal-insulator transition at the criti- 
cal value of the conductance gx, such that samples with 
conductances gx > gr ^'i'^ metals and their conductivity 
remains finite at T ^ while samples with gx < g^ ^-re 
insulators and their conductivity vanishes at T — > 0. 

The main results of our work are as follows: (i) We 
estimate the critical value g^ of the tunnel conductance 
at which the metal-insulator transition in 31? occurs as 



5^ = (l/67r) HEc/6), 



(1) 



where Ec is the charging energy of an isolated grain, (ii) 
We find the expression for the conductivity of a granular 
metal that includes corrections due to Coulomb inter- 
action and holds for all temperatures as long as these 
corrections are small. The corresponding answer can be 
conveniently written separating the correction due to the 
contribution from the large energy scales e > grS from 
that coming from the low energy scales e < grS. Denot- 
ing corrections as Sai and Sa2 respectively we have 



(2a) 



where ctq — '^e^gra^ '^, with a being the size of the single 
grain is the classical Drude conductivity for a granular 
metal (spin included). Correction Sai in Eq. (|2a|l con- 
tains the dimensionality of the array d only as a coeffi- 
cient and is given by the following expression 



6ai 



1 



2ndgT 



In 



grEc 



max (T, grS) 



(2b) 



On the contrary the correction 5a2 in Eq. H2a() that is im- 
portant only at temperatures T < Sgr strongly depends 
on the dimensionality of the array 



Sa2 

CTO 




D^3, 
D = 2, 
D = l. 



(2c) 



Here a = dxx-^/^[l - coth(x) + x/sinh^(x)] « 1.83 
and P = dx x^-^/^ [coth(x) — x/ sinh^(a;)] ~ 3.13 are 
the numerical constants. For a 3D granular system a 
temperature independent term of the order 1 /gT has been 
subtracted in the first line in Eq. if^ . 

Corrections Sai and 6a2 are of a different origin: the 
correction Sai comes from the large energy scales, e > 



gxS where the granular structure of the array dominates 
the physics. The fact that this correction is essentially 
independent of the dimentionality d means that the tun- 
neling of electrons with energies e > gxS can be consid- 
ered as incoherent. On the other hand, correction Sa2 
in Eq. is similar to that obtained for homogeneous 
metals long ago i3] and comes from the low energy scales, 
£ < grS, where the coherent electron motion on the scales 
larger than the grain size a dominates the physics. 

It is important to note that in the low temperature 
regime all temperate dependence of conductivity comes 
from the correction Sa2 ■ At the same time, in this regime 
the correction Sai, though being temperature indepen- 
dent, still exists and can be even larger than Sa2. 

When deriving Eqs. (O we neglected possible weak lo- 
calization corrections that may originate from quantum 
interference of electron waves. This approximation is le- 
gitimate if a magnetic field is applied as in Ref. Ij or 
dephasing is strong due to inelastic processes. 

Now we turn to the description of our model and the 
derivation of Eqs. We consider a c?— dimensional ar- 
ray of metallic grains with the Coulomb interaction be- 
tween electrons. The motion of electrons inside the grains 
is diffusive and they can tunnel from grain to grain. In 
principle, the grains can be clean such that electrons scat- 
ter mainly on grain surfaces. We assume that the sam- 
ple in the absence of the Coulomb interaction would be 
a good metal. For large tunneling conductance we may 
also neglect the nonperturbative charging effects (discret- 
ness of the electron charge) Q, which give an exponen- 
tially small (as exp(— #3^)) contribution to the conduc- 
tivity. Although we assume that the diniensionless tun- 
neling conductance gr is large, it should be still smaller 
than the grain conductance, go, such that gx < go- This 
inequality means that the granular structure is still im- 
portant and the main contribution to the macroscopic 
resistivity comes from the contacts between the grains. 
The limit of the homogeneous metal is reached when the 
grain conductance, go becomes of the order of tunneling 
conductance, 170 ^ gr- 

The system of weakly coupled metallic grains can be 
described by the Hamiltonian 

H = Ha + H, + Y, Uj [^jHn)i'{rj)+i^Hr,)^j{n)], (3) 

where tij is the tunneling matrix element corresponding 
to the points of contact and Vj of i-th and j— th grains. 
The Hamiltonian Hq in Eq. (PJ describes noninteracting 
isolated disordered grains. The term He describes the 
Coulomb interaction inside and between the grains. It 
has the following form 



(4) 



where Cy is the capacitance matrix and fii is the opera- 
tor of electrons number in the i-th grain. In the regime 
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FIG. 1: Diagrams describing the conductivity of granular 
metals: the diagram (a) corresponds to ctq in Eq. Il2all and 
it is the analog of Drude conductivity. Diagrams (b)-(e) de- 
scribing first order correction to the conductivity of granular 
metals due to electron-electron interaction. The solid lines 
denote the propagator of electrons and the dashed lines de- 
scribe effective screened electron-electron propagator. The 
tunneling vertices are described by the circles. The sum of 
the diagrams (b) and (c) results in the conductivity correc- 
tion (5(71 in Eq. 12all . The other two diagrams, (d) and (e) 
result in the correction 5a2- 



under consideration one can neglect the coordinate de- 
pendence of a single grain diffusion propagator. The elec- 
tron hopping between the grains can be included using 
the diagrammatic technique developed in Refs.0,H. 

The conductivity of the granular metals is given by the 
analytical continuation of the Matsubara current-current 
correlator. In the absence of the electron-electron inter- 
action the conductivity is represented by the diagram (a) 
in Fig. 1 that results in high temperature (Drude) con- 
ductivity ao which is defined below Eq. H2a|l . First order 
interaction corrections to the conductivity are given by 
the diagrams (b-e) in Fig. 1. These diagrams are anal- 
ogous to ones considered in Ref. for the correction to 
the conductivity of homogeneous metals. We consider the 
contributions from diagrams (b,c) and (d,e) separately: 
The sum of the diagrams (b,c) results in the following 
correction to the conductivity 



ao 2'KdgT ^ J 



(5) 



where 7 (w) = ^wcoth^, Sq — 257- ( 1 — cos qa) with 
a being the lattice vectors and 
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iuj) (4eq£'c(q) ~ iuj)' 



(6) 



Here the charging energy -Ec(q) — e^/2C(q) is expressed 
in terms to the Fourier transform of the capacitance ma- 
trix C(q). Performing the integration over the frequency 
and summing over the quasimomentum q in Eq. Q with 
the logarithmic accuracy we obtain the correction l|2bp . 
One can see from Eq. © that the contribution Sai in 
Eq. (|2b|) comes from the large energy scales, e > grS 
such that at low temperatures the logarithm is cut off on 
the energy scale grS. 

To obtain the total correction to the conductivity of 
granular metal the two other diagrams, (d) and (e) in 



Fig. 1 should be taken into account. These diagrams 
result in the following contribution to the conductivity 
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(7) 

In contrast to the contribution 5ai in Eq. the main 
contribution to the sum over the quasimomentum q in 
Eq. (O comes from the low momenta, q <^ 1/a. In this 
regime the capacitance matrix, C(q) in Eqs. © and Q) 
has the following asymptotic form 



Substituting Eq. © into Eq. ©, integrating over the 
frequency and summing over the quasimomentum q in 
Eq. (O we obtain the result for the correction 802 in 
Eq. if^ . This type of corrections to the conductivity 
cannot be obtained from the effective AES action. 

Comparing our results in Eqs. (0) with those obtained 
in Ref. '3' using the AES functional we see that the correc- 
tion to the conductivity obtained in Ref. is equivalent 
to the correction 5a i in Eq. (|2a|l . which corresponds in 
our approach to the sum of diagrams (b) and (c) in Fig. 1. 
The correction 5a2 in Eq. (|2a|l becomes important only at 
low temperatures, T < qtS where AES functional is not 
applicable. While in our approach both corrections to the 
conductivity must be small 6ai , Sa2 ^ cto the method of 
Ref. gives a possibility to show that for T 3> grS the 
dependence of the conductivity is logarithmic so long as 

It follows from Eq. that at low temperatures, 

T < gx^, for a 3-D granular array, there are no essen- 
tial corrections to the conductivity coming from the low 
energies since the correction Sa2 is always small. This 
means that the result for the renormalized conductance, 
gx of Ref. for 3D samples within the logarithmic accu- 
racy can be written in the following form 



9t{T) = gr - 7^ In 

DTT 



QtEc 



max {grS, T) 



(9) 



such that it is valid for all temperatures as long as 
the renormalized conductance, gT S> 1. One can 
see from Eq. © that for bare conductance, gr 3> 
{l/&'K)\i\{gTEc / 5) the renormalized conductance, gr is 
always large and the system remains metallic down 
to zero temperatures. In the opposite limit gr < 
(1/ 6tt) ln{gTEc' / S) , the system flows when decreasing the 
temperature to the strong coupling regime, ~ 1 that 
indicates the onset of the insulating phase. We see that 
with the logarithmic accuracy the critical value of the 
conductance g^ is given by Eq. (Q). 

In a similar way we can obtain interaction corrections 
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to the density of states (DOS) of granular metal 



Sryje) 



An ^ 
q 



Im / doj 



tanh[(e - w)/2T] 



iujlAEc{c\)\ 
(10) 

Here vq is the DOS for noninteracting electrons, £q and 
£'c(q) were defined below Eqs. (0) and © respectively. 
Using Eq. (|10|) for a 3D granular array we obtain 
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2ngT 
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Ec9t 



max (T, grb) 



(11a) 



where A 



gra- 



J d^q/ {27:) 



3 p-1 

q • 



For temperatures, 



T ^ qtS the correction to the DOS (|lla|l coincides with 
the one obtained in Ref. 3] using AES approach. It fol- 
lows from Eq. Hlla() that for a 3D array of grains, as in 
case with conductivity, the main contribution to the DOS 
comes from the large energy scales, e > grS- 

Using Eq. H1U|I for a 2D array we obtain the following 
result for the correction to the DOS 



6V2 
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IG^ttt^ 
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In 91^ In ^ 



2 In^ ^ 



T > grS, 
T < gT5. 

(lib) 

Using the relation between the tunneling conductance 
and the diffusion coefhcient, D = gxa^S one can see that 
the temperature dependence of the DOS for T -C grS 
given by Eq. Ijllb|) coincides up to the temperature in- 
dependent constant with the result for the correction to 
the DOS of the homogeneous metal 

A logarithmic dependence of resistivity on tem- 
perature was recently found in high-Tc compounds 
La2-ySryCu04 and BioSr^xLaxCuOe+s hi a very 
strong magnetic field j^, ^ possible granular- 

ity of these samples was suggested in Ref. Re- 
cently the microscopic granularity was directly exper- 
imentally observed in the superconducting state of 
Bi2Sr2CaCu208+s by the STM probe If we accept 
that samples studied in 0, |l3| are indeed microscopically 
granular, we can compare the results of the experiments 
with our predictions. When doing so it is convenient to 
scale three dimensional conductivity to the conductivity 
of CuO planes, a plane- According to our predictions 



dupi 



d\nT 



ane ^ 



(12) 



where the coefficient k — l/27r in the low temperature- 
and k — 1/d in the high temperature regimes. While in 
the low temperature regime the application of Eq. 1)12(1 
is legitimate only under the assumption that electrons 
in different CuO plane are incoherent, in the high tem- 
perature regime the behavior of conductivity according 
to Eq. ((2b|l is logarithmic for any dimension. In this 
regime the real dimensionality d should be replaced by 
d = Z/2, where Z is the (average) number of the contacts 
of each grain with all the adjacent grains. Describing the 



data shown in Fig. 3 of Ref. by our log dependen- 
cies at temperatue T « 5K we extract k ~ 0.4, for Sr 
concentration of y — 0.08 for La2-ySryCu04^ [l3|; for 
the Bi2Sr2-xLaxCuOG-f-s compound we find k ~ 0.2 for 
X = 0.84 La concentration, and k ~ 0.3 for x = 0.76. For 
each particular curve the values k extracted from Fig. 3 
of Ref. [l3| increase with temperature (especially in case 
of LSCO), this is in a complete agreement with our re- 
sults provided that the "coherent-incoherent" crossover 
occurs at about T ~ bK. At higher temperatures k no- 
ticeably exceeds l/27r, supporting the idea of a granular- 
ity of doped cuprates. The questions whether the high 
values of in the superconductors are related to the 
granular structure and what is the origin of this granu- 
larity are beyond the scope of the present work. 

In conclusion, we have investigated transport proper- 
ties of granular metals at large tunneling conductance 
and obtained corrections to the conductivity, Eqs. (0) 
and DOS, Eqs. 1(11(1 due to electron-electron interaction. 
We have shown that at temperatures, T > gxS the gran- 
ular structure of the array dominates the physics. On the 
contrary at temperatures, T < gxS the large-scale coher- 
ent electron motion is crucial. Comparison our results 
with experimental data supports the assumption about 
a granular structure of doped high-Tc cuprates. 
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